The role of multiple marks in epigenetic silencing and the emergence of a stable bivalent chromatin state

We introduce and analyze a minimal model of epigenetic silencing in budding yeast, built upon known biomolecular interactions in the system. Doing so, we identify the epigenetic marks essential for the bistability of epigenetic states. The model explicitly incorporates two key chromatin marks, namely H4K16 acetylation and H3K79 methylation, and explores whether the presence of multiple marks lead to a qualitatively different systems behavior. We find that having both modifications is important for the robustness of epigenetic silencing. Besides the silenced and transcriptionally active fate of chromatin, our model leads to a novel state with bivalent (i.e., both active and silencing) marks under certain perturbations (knock-out mutations, inhibition or enhancement of enzymatic activity). The bivalent state appears under several perturbations and is shown to result in patchy silencing. We also show that the titration effect, owing to a limited supply of silencing proteins, can result in counter-intuitive responses. The design principles of the silencing system is systematically investigated and disparate experimental observations are assessed within a single theoretical framework. Specifically, we discuss the behavior of Sir protein recruitment, spreading and stability of silenced regions in commonly-studied mutants (e.g., sas2, dot1) illuminating the controversial role of Dot1 in the systems biology of yeast silencing.Comment: Supplementary Material, 14 page

Effect of promoter architecture on the cell-to-cell variability in gene expression

According to recent experimental evidence, the architecture of a promoter, defined as the number, strength and regulatory role of the operators that control the promoter, plays a major role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effort that addresses the question of how changes in promoter architecture affect noise in gene expression in a systematic rather than case-by-case fashion. In this article, we make such a systematic investigation, based on a simple microscopic model of gene regulation that incorporates stochastic effects. In particular, we show how operator strength and operator multiplicity affect this variability. We examine different modes of transcription factor binding to complex promoters (cooperative, independent, simultaneous) and how each of these affects the level of variability in transcription product from cell-to-cell. We propose that direct comparison between in vivo single-cell experiments and theoretical predictions for the moments of the probability distribution of mRNA number per cell can discriminate between different kinetic models of gene regulation.Comment: 35 pages, 6 figures, Submitte

Sequential induction of three recombination directionality factors directs assembly of tripartite integrative and conjugative elements

Tripartite integrative and conjugative elements (ICE3) are a novel form of ICE that exist as three separate DNA regions integrated within the genomes of Mesorhizobium spp. Prior to conjugative transfer the three ICE3 regions of M. ciceri WSM1271 ICEMcSym1271 combine and excise to form a single circular element. This assembly requires three coordinated recombination events involving three site-specific recombinases IntS, IntG and IntM. Here, we demonstrate that three excisionases–or recombination directionality factors—RdfS, RdfG and RdfM are required for ICE3 excision. Transcriptome sequencing revealed that expression of ICE3 transfer and conjugation genes was induced by quorum sensing. Quorum sensing activated expression of rdfS, and in turn RdfS stimulated transcription of both rdfG and rdfM. Therefore, RdfS acts as a “master controller” of ICE3 assembly and excision. The dependence of all three excisive reactions on RdfS ensures that ICE3 excision occurs via a stepwise sequence of recombination events that avoids splitting the chromosome into a non-viable configuration. These discoveries expose a surprisingly simple control system guiding molecular assembly of these novel and complex mobile genetic elements and highlight the diverse and critical functions of excisionase proteins in control of horizontal gene transfer

Transcriptional interaction-assisted identification of dynamic nucleosome positioning

<p>Abstract</p> <p>Background</p> <p>Nucleosomes regulate DNA accessibility and therefore play a central role in transcription control. Computational methods have been developed to predict static nucleosome positions from DNA sequences, but nucleosomes are dynamic in vivo.</p> <p>Results</p> <p>Motivated by our observation that transcriptional interaction is discriminative information for nucleosome occupancy, we developed a novel computational approach to identify dynamic nucleosome positions at promoters by combining transcriptional interaction and genomic sequence information. Our approach successfully identified experimentally determined nucleosome positioning dynamics available in three cellular conditions, and significantly improved the prediction accuracy which is based on sequence information alone. We then applied our approach to various cellular conditions and established a comprehensive landscape of dynamic nucleosome positioning in yeast.</p> <p>Conclusion</p> <p>Analysis of this landscape revealed that the majority of nucleosome positions are maintained during most conditions. However, nucleosome occupancy at most promoters fluctuates with the corresponding gene expression level and is reduced specifically at the phase of peak expression. Further investigation into properties of nucleosome occupancy identified two gene groups associated with distinct modes of nucleosome modulation. Our results suggest that both the intrinsic sequence and regulatory proteins modulate nucleosomes in an altered manner.</p

Recombinant Mouse PAP Has pH-Dependent Ectonucleotidase Activity and Acts through A1-Adenosine Receptors to Mediate Antinociception

Prostatic acid phosphatase (PAP) is expressed in nociceptive neurons and functions as an ectonucleotidase. When injected intraspinally, the secretory isoforms of human and bovine PAP protein have potent and long-lasting antinociceptive effects that are dependent on A1-adenosine receptor (A1R) activation. In this study, we purified the secretory isoform of mouse (m)PAP using the baculovirus expression system to determine if recombinant mPAP also had antinociceptive properties. We found that mPAP dephosphorylated AMP, and to a much lesser extent, ADP at neutral pH (pH 7.0). In contrast, mPAP dephosphorylated all purine nucleotides (AMP, ADP, ATP) at an acidic pH (pH 5.6). The transmembrane isoform of mPAP had similar pH-dependent ectonucleotidase activity. A single intraspinal injection of mPAP protein had long-lasting (three day) antinociceptive properties, including antihyperalgesic and antiallodynic effects in the Complete Freund's Adjuvant (CFA) inflammatory pain model. These antinociceptive effects were transiently blocked by the A1R antagonist 8-cyclopentyl-1, 3-dipropylxanthine (CPX), suggesting mPAP dephosphorylates nucleotides to adenosine to mediate antinociception just like human and bovine PAP. Our studies indicate that PAP has species-conserved antinociceptive effects and has pH-dependent ectonucleotidase activity. The ability to metabolize nucleotides in a pH-dependent manner could be relevant to conditions like inflammation where tissue acidosis and nucleotide release occur. Lastly, our studies demonstrate that recombinant PAP protein can be used to treat chronic pain in animal models

Plasmid pP62BP1 isolated from an Arctic Psychrobacter sp. strain carries two highly homologous type II restriction-modification systems and a putative organic sulfate metabolism operon

The complete nucleotide sequence of plasmid pP62BP1 (34,467 bp), isolated from Arctic Psychrobacter sp. DAB_AL62B, was determined and annotated. The conserved plasmid backbone is composed of several genetic modules, including a replication system (REP) with similarities to the REP region of the iteron-containing plasmid pPS10 of Pseudomonas syringae. The additional genetic load of pP62BP1 includes two highly related type II restriction-modification systems and a set of genes (slfRCHSL) encoding enzymes engaged in the metabolism of organic sulfates, plus a putative transcriptional regulator (SlfR) of the AraC family. The pP62BP1 slflocus has a compact and unique structure. It is predicted that the enzymes SlfC, SlfH, SlfS and SlfL carry out a chain of reactions leading to the transformation of alkyl sulfates into acyl-CoA, with dodecyl sulfate (SDS) as a possible starting substrate. Comparative analysis of the nucleotide sequences of pP62BP1 and other Psychrobacter spp. plasmids revealed their structural diversity. However, the presence of a few highly conserved DNA segments in pP62BP1, plasmid 1 of P. cryohalolentis K5 and pRWF-101 of Psychrobacter sp. PRwf-1 is indicative of recombinational shuffling of genetic information, and is evidence of lateral gene transfer in the Arctic environment

The Waddlia Genome: A Window into Chlamydial Biology

Growing evidence suggests that a novel member of the Chlamydiales order, Waddlia chondrophila, is a potential agent of miscarriage in humans and abortion in ruminants. Due to the lack of genetic tools to manipulate chlamydia, genomic analysis is proving to be the most incisive tool in stimulating investigations into the biology of these obligate intracellular bacteria. 454/Roche and Solexa/Illumina technologies were thus used to sequence and assemble de novo the full genome of the first representative of the Waddliaceae family, W. chondrophila. The bacteria possesses a 2′116′312bp chromosome and a 15′593 bp low-copy number plasmid that might integrate into the bacterial chromosome. The Waddlia genome displays numerous repeated sequences indicating different genome dynamics from classical chlamydia which almost completely lack repetitive elements. Moreover, W. chondrophila exhibits many virulence factors also present in classical chlamydia, including a functional type III secretion system, but also a large complement of specific factors for resistance to host or environmental stresses. Large families of outer membrane proteins were identified indicating that these highly immunogenic proteins are not Chlamydiaceae specific and might have been present in their last common ancestor. Enhanced metabolic capability for the synthesis of nucleotides, amino acids, lipids and other co-factors suggests that the common ancestor of the modern Chlamydiales may have been less dependent on their eukaryotic host. The fine-detailed analysis of biosynthetic pathways brings us closer to possibly developing a synthetic medium to grow W. chondrophila, a critical step in the development of genetic tools. As a whole, the availability of the W. chondrophila genome opens new possibilities in Chlamydiales research, providing new insights into the evolution of members of the order Chlamydiales and the biology of the Waddliaceae

Intronic L1 Retrotransposons and Nested Genes Cause Transcriptional Interference by Inducing Intron Retention, Exonization and Cryptic Polyadenylation

Transcriptional interference has been recently recognized as an unexpectedly complex and mostly negative regulation of genes. Despite a relatively few studies that emerged in recent years, it has been demonstrated that a readthrough transcription derived from one gene can influence the transcription of another overlapping or nested gene. However, the molecular effects resulting from this interaction are largely unknown.Using in silico chromosome walking, we searched for prematurely terminated transcripts bearing signatures of intron retention or exonization of intronic sequence at their 3' ends upstream to human L1 retrotransposons, protein-coding and noncoding nested genes. We demonstrate that transcriptional interference induced by intronic L1s (or other repeated DNAs) and nested genes could be characterized by intron retention, forced exonization and cryptic polyadenylation. These molecular effects were revealed from the analysis of endogenous transcripts derived from different cell lines and tissues and confirmed by the expression of three minigenes in cell culture. While intron retention and exonization were comparably observed in introns upstream to L1s, forced exonization was preferentially detected in nested genes. Transcriptional interference induced by L1 or nested genes was dependent on the presence or absence of cryptic splice sites, affected the inclusion or exclusion of the upstream exon and the use of cryptic polyadenylation signals.Our results suggest that transcriptional interference induced by intronic L1s and nested genes could influence the transcription of the large number of genes in normal as well as in tumor tissues. Therefore, this type of interference could have a major impact on the regulation of the host gene expression

Analysis of the Neurotoxin Complex Genes in Clostridium botulinum A1-A4 and B1 Strains: BoNT/A3, /Ba4 and /B1 Clusters Are Located within Plasmids

BACKGROUND: Clostridium botulinum and related clostridial species express extremely potent neurotoxins known as botulinum neurotoxins (BoNTs) that cause long-lasting, potentially fatal intoxications in humans and other mammals. The amino acid variation within the BoNT is used to categorize the species into seven immunologically distinct BoNT serotypes (A-G) which are further divided into subtypes. The BoNTs are located within two generally conserved gene arrangements known as botulinum progenitor complexes which encode toxin-associated proteins involved in toxin stability and expression. METHODOLOGY/PRINCIPAL FINDINGS: Because serotype A and B strains are responsible for the vast majority of human botulism cases worldwide, the location, arrangement and sequences of genes from eight different toxin complexes representing four different BoNT/A subtypes (BoNT/A1-Ba4) and one BoNT/B1 strain were examined. The bivalent Ba4 strain contained both the BoNT/A4 and BoNT/bvB toxin clusters. The arrangements of the BoNT/A3 and BoNT/A4 subtypes differed from the BoNT/A1 strains and were similar to those of BoNT/A2. However, unlike the BoNT/A2 subtype, the toxin complex genes of BoNT/A3 and BoNT/A4 were found within large plasmids and not within the chromosome. In the Ba4 strain, both BoNT toxin clusters (A4 and bivalent B) were located within the same 270 kb plasmid, separated by 97 kb. Complete genomic sequencing of the BoNT/B1 strain also revealed that its toxin complex genes were located within a 149 kb plasmid and the BoNT/A3 complex is within a 267 kb plasmid. CONCLUSIONS/SIGNIFICANCE: Despite their size differences and the BoNT genes they contain, the three plasmids containing these toxin cluster genes share significant sequence identity. The presence of partial insertion sequence (IS) elements, evidence of recombination/gene duplication events, and the discovery of the BoNT/A3, BoNT/Ba4 and BoNT/B1 toxin complex genes within plasmids illustrate the different mechanisms by which these genes move among diverse genetic backgrounds of C. botulinum

Run-Off Replication of Host-Adaptability Genes Is Associated with Gene Transfer Agents in the Genome of Mouse-Infecting Bartonella grahamii

The genus Bartonella comprises facultative intracellular bacteria adapted to mammals, including previously recognized and emerging human pathogens. We report the 2,341,328 bp genome sequence of Bartonella grahamii, one of the most prevalent Bartonella species in wild rodents. Comparative genomics revealed that rodent-associated Bartonella species have higher copy numbers of genes for putative host-adaptability factors than the related human-specific pathogens. Many of these gene clusters are located in a highly dynamic region of 461 kb. Using hybridization to a microarray designed for the B. grahamii genome, we observed a massive, putatively phage-derived run-off replication of this region. We also identified a novel gene transfer agent, which packages the bacterial genome, with an over-representation of the amplified DNA, in 14 kb pieces. This is the first observation associating the products of run-off replication with a gene transfer agent. Because of the high concentration of gene clusters for host-adaptation proteins in the amplified region, and since the genes encoding the gene transfer agent and the phage origin are well conserved in Bartonella, we hypothesize that these systems are driven by selection. We propose that the coupling of run-off replication with gene transfer agents promotes diversification and rapid spread of host-adaptability factors, facilitating host shifts in Bartonella